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Abstract. We present the results of a global neutrino oscillation data analysis 
within the three-flavour framework. We include latest results from the MINOS long- 
fS |. baseline experiment (including electron neutrino appearance as well as anti-neutrino 
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data), updating all relevant solar (SK II+III), atmospheric (SK I+II+III) and reactor 
. (KamLAND) data. Furthermore, we include a recent re-calculation of the anti-neutrino 

i-Sh I fluxes emitted from nuclear reactors. These results have important consequences for the 

analysis of reactor experiments and in particular for the status of the mixing angle #13. 
, In our recommended default analysis we find from the global fit that the hint for non- 

zero 0i3 remains weak, at 1.8cr for both neutrino mass hierarchy schemes. However, we 
■ discuss in detail the dependence of these results on assumptions concerning the reactor 

t N *»- \ neutrino analysis. 
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d ■ 1. Introduction 

The discovery of neutrino mixing and oscillations provides firm evidence for physics 
beyond Standard Model, opening a new era in particle physics. Here we update the 
three-neutrino oscillation results of Ref. |l| |i| with a special emphasis on the new 
reactor anti-neutrino flux results of Refs. p, |6( which have an important impact on 
the determination of the mixing angle 6*13. We include new data from the MINOS 
Collaboration, both for — > z/ f transitions 0] and disappearance 0] , the latest 
Super-Kamiokande (SK) solar IOL LL1[ and atmospheric 12| neutrino data, as well as 



recent KamLAND reactor data [13j. Our goal is to summarize the results of the three- 
flavour neutrino oscillation analysis paying attention to sub-leading three-flavor effects 
where they are most relevant, as well as to the effects of the new anti-neutrino flux 

% When referring to this paper we include also its updates available at arXiv:0808.2016v3 [hep-ph]. 
Further technical details of the analysis as well as earlier experimental references are also given in our 
previous review in Q • For other global analyses see 0, 0] . 
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emitted from nuclear reactors which has been re-evaluated in Ref. [5|. The reported 
value of the v e flux is about 3% higher than from previous calculations. This has 
important consequences for the interpretation of data from reactor experiments. We 
discuss the implications of the new reactor neutrino flux for the determinations of 
oscillation parameters, in particular its effect on the mixing angle #13. We find that 
due to the new fluxes the results depend on the inclusion of short-baseline reactor data 
from distances < 100 m. 

In Sec. [2]we present the updated analysis in the "atmospheric sector" , discussing the 



Super Kamiokande I+II+III data [12| in Sec 12.11 and the MINOS disappearance results 
taking into account the solar squared-mass splitting and discussing the slight tension 
between neutrino and anti-neutrino data in Sec. 12.21 In Sec. I2.3l we focus on the combined 
SK+MINOS analysis and the determination of #13 from these data. Sec. |3] contains the 
discussion of the reactor neutrino data in the light of the new predicted anti-neutrino 
fluxes, and in Sec. |4]we discuss solar neutrino data as well as KamLAND and the other 
reactor experiments. The results of the global fit are summarized in Sec. [5] including a 
detailed discussion of the status of 6*13. In particular we discuss the dependence of the 
#13 determination upon assumptions concerning the reactor anti-neutrino data analysis. 
Conclusions follow in Sec. |6j 

2. The atmospheric sector - SK I+II+III and MINOS 

2.1. Super- Kamiokande I+II+III data 

We include in our analysis the full sample of atmospheric neutrino data from all three 
phases of the Super-Kamiokande experiment 12|, using directly the x 2 map provided 



by the Super-Kamiokande collaboration. The atmospheric neutrino oscillation analysis 
is performed within the one mass scale approximation, neglecting the effect of the solar 
mass splitting, as in our previous papers jl|, 0] @- 

In Fig. [1] we show the results of our previous atmospheric neutrino analysis from 
Refs. [l|, Q] including data from the first phase of the Super-Kamiokande experiment, 
compared to the new analysis included in this update, based on the Super-Kamiokande 
data from all three phases of the experiment. The differences between the two 
analyzes are clearly noticeable, with an improved determination of both oscillation 
parameters. As we will see in the following, once data from neutrino disappearance 
at the MINOS long-baseline experiments are included in the analysis, the improvement 
in the determination of Am^ due to the recent atmospheric neutrino data is "hidden" 
by the more constraining restrictions imposed by long-baseline data. Nevertheless, this 
new atmospheric analysis will be important when constraining the mixing angle #23 and 
also #13. 



§ Preliminary results towards a full 3-flavour atmospheric neutrino analysis have been presented in [14 
we look forward to the corresponding information becoming publicly available. 
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Figure 1. Comparison of our previous SK-I atmospheric analysis (in black) and the 
new analysis using the latest atmospheric data from SK-I, SK-II, SK-II (in red) for 
inverted mass hierarchy. Best fit points denoted by a dot or star follow the same colour 
code. 



> 
CD 

CO 

O 



CVJ CO 



r~i — r~i — r~r~i — i i | i i — r~r~| — r~i — r~r 



MINOS v 3.4x10" p.o.t. 
MINOS v 7.2xl0 2 ° p.o.t. 

11 r 




_L 



0.25 



0.5 

sin 9 



0.75 



23 




1 



MINOS v only 
MINOS v^only 
MINOS v + v 



Figure 2. Determination of atmospheric neutrino parameters by the MINOS long- 
baseline experiment. Left panel: previous (3.36 x 10 20 p.o.t.) versus current (7.2xl0 20 
p.o.t.) data. Right panel: allowed regions from recent MINOS (7.2xl0 20 p.o.t.) 
data, using neutrinos-only, anti-neutrinos-only, and the combination. Both panels 
assume normal hierarchy. As before, best fit points follow the same colour code. 



2.2. MINOS disappearance data 

At the Neutrino 2010 Conference, the MINOS Collaboration has presented new data 
from their searches of disappearance, both from the neutrino (7.2 xlO 20 p.o.t.) and 
the anti-neutrino (1.71 xlO 20 p.o.t) running mode jsj. We perform a re-analysis of the 
data from js| within a full three-flavour framework using the GLoBES software 15 



to simulate the experiment. In addition to matter effects, we include also the effect 
of A well as 8± 3 and the CP-phase 5 in the analysis of the disappearance and 

appearance channels. Since in our analysis MINOS data are the only ones sensitive to 
the phase 5 we always minimize the MINOS x 2 with respect to 5. 
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In Fig.[2]we compare the analysis of new MINOS data with the previous data release 
(left panel). Apart from the smaller size of the allowed regions due to the increase in 
statistics, we notice also a change in the shape of the regions. This follows from the 
inclusion of the sub-leading three-flavour effects included in this new analysis. 

In the right panel we illustrate the impact of the (still small) anti-neutrino sample. 
From the figure one can see that there is a slight tension between the neutrino and 
anti-neutrino results: there is no overlap of the allowed regions at less than 90% CL. 
However, at 3a the results of both are fully consistent. We find the following \ 2 minima 
and goodness-of-fit (GOF) values: 

Xmin,, = 24.4/(27 - 2) GOF = 49.6% 
v: ^ = 15.0/(13-2) GOF = 18.4% (1) 
v + v-- xLn,tot = 46.1/(40 - 2) GOF = 17.3% 
Hence the combined neutrino and anti-neutrino fit provides still an acceptable GOF. 
Using the consistency test from Ref. (l6| yields Xpg = Xmin.tot ~ Xmi™ - Xmm,p = 6 - 6 - 



The value of Xpg m ust be evaluated for 2 degrees of freedom, which implies that neutrino 
and anti-neutrino data are consistent with a probability of 3.7%. This number indicates 
a slight tension between the sets, at the level of about 2.1<r. In the following we will use 
only neutrino data in the global analysis. It is clear from Fig. [2] (right) that adding also 
anti- neutrino data would have negligible impact on the global result. 

2.3. The atmospheric sector: combined MINOS + atmospheric analysis and #13 

Combining the new atmospheric and MINOS disappearance data we obtain new global 
constraints on the atmospheric neutrino oscillation parameters. The results are shown in 
Fig-El As before, Am^ is determined mainly by the MINOS data, while the atmospheric 
data are more important in determining the mixing parameter sin 2 #23 • Notice that the 
sub- leading effects of Am^ lead to different best fit points for jAm^l, that depend on 
whether the mass hierarchy is normal (NH) or inverted (IH). We find the following best 
fit values with errors at la: 

I Am* I- J 2 - 45±0 - 09 xKTW (NH) 

1 311 ~ [ 2.34+°;J° xlO- 3 eV 2 (IH) [Z) 

The corresponding allowed regions are shown in Fig. [3] (right). The reason for this 
apparent shift is just a result of our parameterization, using Am^ for both hierarchies, 
and changing only the sign of it to distinguish them. Hence, in NH the "largest" 
frequency is given by lAmjjjJ, while in IH it is lAmU + Am 2 1; which explains why 
(Am^l is smaller for IH. Eq. [2] shows that these sub- leading effects must be included 
given the present accuracy, since they are at the level of the la error on Am^. 

Now we move to appearance MINOS data. The MINOS Collaboration has recently 
also reported new data from the search of — > v e transitions in the Fermilab NuMI 



)20 

more than twice the size of the previous data release [17|. The new MINOS far detector 



beam [7J. The new data are based on a total exposure of 7x10 protons-on-target, 
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Figure 3. Determination of the atmospheric oscillation parameters. Left: interplay 
of atmospheric (black) and MINOS disappearance (blue) data and the combination 
(red/shaded region) for normal hierarchy at 90% CL (dashed) and 3cr (solid). Right: 
combined allowed regions for normal (black curves) and inverted hierarchy (colored 
regions) at 90%, 95%, 99%, 99.73% CL. 



data consists of 54 electron neutrino events, while, according to the measurements in the 
MINOS Near Detector, 49.1 ± 7.0(stat) ± 2.7 (syst) background events were expected. 
Hence the observed number of events is in agreement with background expectations 
within 0.7o" and the hint for a non-zero value of #13 present in previous data [13] has 
largely disappeared. In fact, we see that once we include the new MINOS data in our 
analysis, a smaller best fit point of 6* 13 is obtained and, as a result, the hint for #13 is 
less significant than before: for both hierarchies we find only a 0.8a hint when using 
new MINOS data versus 1.3a obtained with the previous MINOS appearance data, see 



e.g. 



18 for a discussion. 



Atmospheric neutrino data from Super-Kamiokande I+II+III described in the 



previous section implies a best fit point very close to 9 i3 = [12J, with Ax = 0.0(0.3) 
for #13 = for NH (IH). However, in the combination with MINOS disappearance and 
appearance data we even find a slight preference for 6>i 3 > 0, with A% 2 = 1.6(1.9) at 
#13 = for NH (IH). As shown in Fig. [4] this happens due to a small mismatch of the 
best fit values for [Am^l at #13 = 0, which can be resolved by allowing for non-zero 
values of #13 j^]. This is similar to the hint for > coming from a slight tension 
between solar and KamLAND data, see Ref. [l|. Therefore, the hint for #13 > from 
atmospheric + LBL data becomes slightly stronger compared to the previous data. 

3. New reactor fluxes and implications for oscillation parameters 

Up to very recently the interpretation of neutrino oscillation searches at nuclear power 



plants was based on the calculations of the reactor v e flux from Ref. [19[. Indeed, the 
observed rates at all reactor experiments performed so-far at distances L < 1 km are 
consistent with these fluxes, therefore setting limits on v e disappearance. Recently the 



Global neutrino data and recent reactor fluxes: status of three-flavour oscillation parameters^ 




Figure 4. Allowed regions at ler and 90% CL for atmospheric (SK) and MINOS 
disappearance and appearance data in the plane of sin 2 #13 and Am 2 ^ for NH (left) and 
IH (right). Combined data is shown as shaded/red region at 90% CL. The black dot, 
blue star, white diamond correspond to the best fit points of SK, MINOS, SK+MINOS, 
respectively. 



flux of v e emitted from nuclear power plants has been re-evaluated [5[ , yielding roughly 
3% higher neutrino fluxes than assumed previously. As discussed in Ref. jfl] this might 
indicate an anomaly in reactor experiments at L < 1 km, which according to the new 
fluxes observe a slight deficit. For the Chooz and Palo Verde experiments at L ~ 1 km 
a non-zero 6*13 could lead to v e disappearance accounting for the reduction of the rate. 
However, Am 2 3 and #13 driven oscillations will have no effect in short-baseline (SBL) 
experiments with L < 100 m. 

Motivated by this situation we include here also the SBL reactor experiments 



Bugey4 |2l|, ROVNO [22 



251 ] . and Gosgen [26 



Bugey3 [23j, Krasnoyarsk [24j, ILL 
via the rate measurements summarized in Table II of p], in addition to the fit of the 
KamLAND, Chooz 27[ and Palo Verde 28] experiments. We use the neutrino fluxes 
from the isotopes 235 U, 239 Pu, 238 U, 241 Pu obtained in j^J. For each reactor experiment 
we take into account the appropriate relative contribution of the isotopes to the total 
flux and we include the uncertainty on the integrated flux for each isotope given in 
Table I of (sf, correlated among all experiments. The total error on fluxes from 235 U, 
239 Pu, 241 Pu are at the level of 2%, where we assume that an error of 1.8% is fully 
correlated among the three isotopes, due to a common normalization uncertainty of the 



corresponding beta-spectra measured in [19 



The SBL reactor data is summarized in Fig. [5j We show the observed rate relative 
to the predicted rate based on old and new flux calculations. Due to the slightly higher 
fluxes according to jsj all experiments observe a smaller ratio with the new fluxes. In 
Fig. [5] we show also the result of a fit to the data with the predicted fluxes, allowing the 
four neutrino fluxes to float in the fit subject to the uncertainties as described above. In 
the fit we assume that the experimental systematic errors of the three data points from 
Bugey3, Gosgen, and Krasnoyarsk, as well as Bugey4 and ROVNO are correlated, due 
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Figure 5. Short-baseline reactor data. We show the observed rate relative to the 
predicted rate based on old [l^] (blue) and new (red) flux calculations. Small error 
bars show statistical and uncorrclated systematic uncertainties, large error bars include 
in addition correlated systematic Al uncertainties. The solid histograms correspond to 
the fitted prediction shifted due to the uncertainty on the fluxes |6| , as indicated in the 
right panel. The dashed lines show the best fit assuming a free overall normalization 
of reactor fluxes. See text for details. 



to the same experimental technique. We obtain x 2 — 8.1(13.0) for 12 degrees of freedom 
using old (new) fluxes. Clearly old fluxes provide a better fit to the data, whereas the 
X 2 for new fluxes is still acceptable (P- value of 37%). Such a good fit can be obtained 
by a rescaling of the fluxes (subject to the quoted uncertainties) as shown in the right 
panel. 

The dashed lines in the figure correspond to a fit where we introduce an overall 
factor / in front of the fluxes, which we let float freely in the fit. For the old fluxes we 
find the best fit value of / = 0.984 with / = 1 within the la range. In contrast, for 
new fluxes we obtain / = 0.942 ± 0.024, and / = 1 disfavored with A% 2 = 6.2 which 
corresponds to about 2.5a. This is the origin of the "reactor anti-neutrino anomaly" 
discussed in jsj] . A possible explanation of this anomaly could be the presence of a sterile 
neutrino "visible" in this oscillation channel but not in the solar and/or atmospheric 
conversions, the so-called 3+1 scenario, see e.g. {29|. However, within the uncertainties 
on the neutrino flux prediction, the goodness of fit of the new fluxes to SBL reactor 
data is still rather good. Given this somewhat ambiguous situation, in the following we 
will present results for 3-flavour oscillations adopting different assumptions on reactor 
neutrino fluxes: (a) motivated by the excellent goodness of fit of SBL data to the new 
flux prediction we take fluxes and the quoted uncertainties at face value, and (b) we 
introduce the free flux normalization / in the fit. This second option takes into account 
the possible presence of a sterile neutrino or some other correlated effect on all reactor 
neutrino fluxes. In this scenario the SBL reactor experiments effectively serve as near 
detectors determining the flux which is then used as input for the oscillation analysis at 
longer baselines. 
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Figure 6. Left: A% 2 as a function of sin 2 #13 for the Chooz and Palo Verde 
(PV) reactor experiments (green), and in combination with the short-baseline (SBL) 
experiments from Fig. [5j Solid (dashed) curves refer to the new (old) reactor anti- 
neutrino fluxes. For the blue curves ("free norm") a free overall normalization factor 
has been introduced for the reactor fluxes. In this figure we fix Am^ = 2.45 x 10 -3 eV 2 . 
Right: Contours in the plane of sin 2 #13 and the flux normalization /. Colored regions 
(curves) correspond to Chooz + PV with (without) including the SBL experiments. 

We show the A% 2 from the Chooz and Palo Verde experiments as a function of 
sin 2 9 13 in Fig. [6] (left) for various assumptions on the fluxes. If the new fluxes are 
taken at face value and SBL reactor experiments are not included in the fit (solid green 
curve), we obtain from Chooz and Palo Verde a hint for # 13 > at about 90% CL, 
with a best fit value at sin 2 6*13 = 0.021. In this case v e disappearance due to #13 
accounts for the suppression of the observed rate at L ~ 1 km relative to the slightly 
increased prediction from the new fluxes. However, as soon as SBL reactor experiments 
are included in the fit, the hint essentially disappears and #13 = is consistent within 
la. This can be understood from Fig. [5j which shows that the SBL reactor experiments 
pull down the flux predictions, leaving less room for a suppression at 1 km due to #13. 
The upper limit on sin 2 # 13 is very similar with old and new fluxes when SBL data are 
included, irrespective of whether the normalization is left free or not. These results are 
in agreement with Ref. (sf. Fig. |6] (right) shows the correlation between sin 2 #13 and the 
flux normalization / with and without SBL experiments. 

To summarize this section we emphasize that due to the tension between the new 
flux predictions and the short baseline oscillation data, including or not the SBL reactors 
in the fit leads to different results concerning the extracted value of #13. This will indeed 
be seen in our subsequent fit results, see for example, Sec. [5j 

4. Solar + KamLAND analysis in the light of new reactor fluxes 



The release of new solar data from the second and third phase of the Super- 
Kamiokande 10|,Qj| experiment, new data from the reactor experiment KamLAND [13| |. 
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Figure 7. Left: 2a allowed regions in the sin 2 612 — Am|i plane from the analysis of 
solar + KamLAND data minimizing over $13. The different curves show the results 
obtained using different assumptions for reactor data, as indicated in the legend. For 
comparison we show as grey-shaded area the region obtained in our previous solar + 
KamLAND data analysis. Right: region allowed at 90%, 95%, 99%, 99.73% CL in our 
recommended analysis of solar + KamLAND data including SBL reactor results. 



and the new predictions for the reactor anti-neutrino fluxes require a full revision of the 
solar + KamLAND neutrino data analysis. In Fig. [7] we show the la allowed region 
from the solar + KamLAND neutrino data using different assumptions in the reactor 
data analysis, such as including or not including the short-baseline data or the use of a 
free normalization factor for the reactor anti-neutrino fluxes. In all cases the shift of the 
allowed region is not very significant, and the best fit point values vary between 0.304 and 
0.320 for sin 2 # 12 while the solar mass splitting Am^i goes from 7.59 to 7.64 x 10~ 5 eV 2 . 
The variation in sin 2 d\i is a bit larger because this parameter is correlated with the 
shift on sin 2 #13, as shown in Fig. El 

In the left panel of Fig. El we show the 2a allowed regions in the sin 2 #1 2 — sin 2 # 13 
plane for different choices of the reactor data analysis. For comparison we also show in 
grey the la allowed region obtained from solar data only. In the right panel of Fig. El 
we show the constraints on #13 from the combination of solar and KamLAND data. 
For the global neutrino oscillation fit presented in the following we will use the analysis 
of KamLAND + SBL data without free normalization, labelled as "sol+KL+SBL" in 
Fig. El In this case we get the following best fit value for 6*13 : 

sin 2 13 = 0.023+°;^ (solar + KamLAND) (3) 

with Ax 2 (sin 2 6 l i3 = 0) = 2.9, and therefore a 1.7<r hint for #13 ^ coming from the 
solar sector. Comparing with our previous analysis the inclusion of new solar and 
KamLAND data and the new reactor fluxes results in a similar best fit value for the (9 13 
mixing angle (before we got sin 2 #13 = 0.022), but a slightly larger significance for non- 
zero 813 (before: A% 2 = 2.2). The origin for this is mainly the preference of KamLAND 
data for a non-zero #13 visible in the left panel of Fig. El For KamLAND #13 acts mainly 




Figure 8. Left: 2cr allowed region in the sin 2 Q\% ~ sm2 ^13 plane obtained from 
KamLAND data using different assumptions on the reactor data analysis. For 
comparison we also show the 2 a allowed region from solar data only (grey-shaded 
area). Right: A\ 2 as a function of sin 2 #13 for the solar + KamLAND data analysis 
under the same assumptions as in the left panel. 



as an overall reduction of the spectrum, and therefore the increased event rate due to 
the new reactor fluxes can be compensated by a non-zero #13 in KamLAND. 

Let us mention that Ref. jfi] also investigates the implication for #13 of reactor 
neutrino data in the light of the new fluxes. Our results are in reasonable agreement, 
although minor quantitative differences occur presumably due to different data used as 
well as different analysis strategies. 



5. Global 3-neutrino analysis and status of 13 

Let us now present the results of the global analysis combining all data mentioned in 
the previous sections [jj As default for the reactor analysis we use the new anti- neutrino 
flux predictions and include in the analysis the SBL reactor experiments as discussed in 
Sec. El 

Fig. M shows the x 2 profile as a function of sin 2 #13 for various data samples. In the 
upper part of Table [1] we display the corresponding best fit values and the significance 
for 9 13 > 0. In our standard recommended analysis (new reactor fluxes, SBL reactors 
included) we find no significant hint from Chooz and Palo Verde data, but the 1.7a hint 
from solar + KamLAND (see Sec. HJ) and the 1.3a (1.4<r for IH) hint from atmospheric 
+ MINOS data (see Sec. 12.3ft combine to a global hint at 1.8a for NH and IH, to 
be compared with the 1.5a obtained in our previous analysis. In the lower part of 
Table CD we discuss how this result depends on details of the reactor neutrino analysis. 
If SBL data are not used in the fit the significance for 6> 13 > is pushed close to 3a 
because as discussed in Sec. [3] in this case also Chooz and Palo Verde prefer 6*13 > at 
about 90% CL. However, in the flux-free reactor analysis as well as with the old reactor 

|| Sec also Refs. [HQ for previous experimental references. 
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Figure 9. Constraint on sin 2 6*13 from different data sets, shown for NH (left) and 
IH (right). The curves labeled "CH+PV+SBL" include the Chooz, Palo Verde and 
the short-baseline reactor experiments, "solar+KL+SBL" include solar, KamLAND 
and short-baseline reactor data, and "atm + LBL" include Super-K atmospheric data, 
MINOS (disappearance and appearance), and K2K. The results from our previous 2010 
analysis are also shown for comparison. 




Figure 10. Illustration of the interplay of the global data on the sin 2 813 bound. Left: 
NH and right: IH. 



fluxes the hint decreases to about 1.4a and 1.8a, respectively. The entry in the table 
labeled "global without reactors" comes from atmospheric and solar neutrinos plus data 
from the MINOS long-baseline experiment. Therefore, these results are independent 
of any ambiguity due to reactor fluxes, and we observe that a non-trivial limit on 6* 13 
emerges even in this case. Let us mention that here we always assume the AGSS09 solar 
model 30 1 . As discussed previously [H, Q there is a minor dependence of the hint for 
9 1 3 on this assumption. 

Fig. [10] illustrates the interplay of the various data sets in the plane of sin 2 6*13 and 
Am^. In Table [2] we summarize the determination of neutrino oscillation parameters 
for our reference default reactor analysis. As discussed in Sec. H]the impact of this choice 
upon the leading oscillation parameters is small. We find that inverted hierarchy gives a 
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Table 1. The best fit values for sin 2 #13 with la errors, the significance of the #13 > 
hint, and the upper bound on sin 2 #13 at 3er for different data samples and for different 
reactor neutrino data assumptions. For a given global analysis the upper (lower) 
numbers refer to normal (inverted) neutrino mass hierarchy. We always use the new 
reactor fluxes 0] except for the row labeled "old fluxes" which uses previous results [l^] . 
The row labeled "free norm" assumes a free reactor anti-neutrino flux normalization. 



parameter 


best fit ±1(7 


2(7 


3(7 


Am|i [10" 5 eV 2 ] 


7 c;q+o.20 

' - OJ -0.18 


7.24-7.99 


7.09-8.19 


Am 2 ! [10~ 3 eV 2 ] 


2.45 ± 0.09 
-(2-341^) 


2.28 - 2.64 
-(2.17-2.54) 


2.18 - 2.73 
-(2.08 - 2.64) 


sin 2 0i2 


n 01 9+0.017 
U.OlZ_ 015 


0.28-0.35 


0.27-0.36 


sin 2 23 


0.51 ±0.06 
0.52 ±0.06 


0.41-0.61 
0.42-0.61 


0.39-0.64 


sin 2 13 


n Q1Q+0.009 
U.U1U_ 006 

1 3+ - 009 
u.u±o_ o 7 


< 0.027 

< 0.031 


< 0.035 

< 0.039 



Table 2. Neutrino oscillation parameters summary. For Am| l5 sin 2 #23, and sin 2 #13 
the upper (lower) row corresponds to normal (inverted) neutrino mass hierarchy. Wc 
assume the new reactor anti-neutrino fluxes [H] and include short-baseline reactor 
neutrino experiments in the fit. 



slightly better fit, however, with only A% 2 = 0.54 with respect to the best fit in normal 
hierarchy. 



Global neutrino data and recent reactor fluxes: status of three-flavour oscillation parameter s!3 



6. Summary and conclusions 

We have presented an updated global fit to world's neutrino oscillation data. The 
recent re-evaluation of the anti-neutrino fluxes emitted in nuclear power plants j^] 
introduces some ambiguity in the results obtained for the mixing angle #13. Since 
the new predictions are in slight disagreement with data from short-baseline (SBL) 
reactor experiments, with L < 100 m, it becomes necessary to include these data in 
the fit. A flux-free analysis of SBL data prefers an off-set of the reactor neutrino flux 
of about 6% from the predicted value with a significance of about 2.5<r. Taken at face 
value this might indicate either some un-accounted systematic effect (either in the new 
calculations or in the reactor data), or even the presence of some kind of new physics such 
as sterile neutrino oscillations with Am 2 ~ 1 eV 2 I6fl. Here we stick to the three-flavour 



framework, the sterile neutrino hypothesis will be discussed elsewhere [311 • 

Despite this hint for a deviation of the observed reactor anti-neutrino flux from its 
prediction, the goodness-of-fit of the SBL reactor neutrino data with the new fluxes is 
still very good (x 2 = 13 for 12 degrees of freedom). Motivated by this result, we adopt 
as our recommended default analysis the new fluxes and include the SBL data in the 
fit. In such a way we obtain a hint for #13 > at 1.8a, coming from a preference for a 
finite #13 from KamLAND data combined with a somewhat weaker hint from the joint 
analysis of atmospheric + MINOS data. In Table Q] we show in detail how the global 
result depends on the assumptions on the reactor neutrino analysis, yielding hints for 
6*i3 > ranging between 1.4<r and 2.8<r, with best fit values between sin 2 6*13 = 0.007 
and 0.027. This somewhat ambiguous situation regarding # 13 emerges from the slight 
tension due to the new reactor flux predictions with data. It will be interesting to see 
how the upcoming results from new reactor as well as accelerator experiments searching 



for # 13 will contribute to resolving the issue, see Ref. [18[ for an overview and references. 

The main results of our recommended default analysis of three-neutrino oscillation 
parameters are summarized in Table. [2] and in the right panels of Figs. |3] and [7] for the 
leading "atmospheric" and "solar" oscillation parameters, as well as in Fig. [10] for the 
mixing angle #13. 
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